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Context. We present stellar yields computed from detailed models of intermediate mass asymptotic giant branch stars of low metallicity. In 
this work, the whole main microphysics inputs have been updated, and in particular a-enhancement is explicitly taken into account both in the 
Q . opacities and equation of state. 

Aims. The target of this work is to provide a basis to test the reliability of the AGB self-enrichment scenario for Globular Clusters of 
C/3 intermediate metallicity. These Globular Clusters exhibit well defined abundance patterns, which have often been interpreted as a consequence 
of the pollution of the interstellar medium by the ejecta of massive AGBs. 

Methods. We calculated a grid of intermediate mass models with metallicity Z=0.001; the evolutionary sequences are followed from the 
, pre-Main sequence along the whole AGB phase. We focus our attention on those elements largely studied in the spectroscopic investigations 

of Globular Clusters stars, i.e. oxygen, sodium, aluminum, magnesium and fluorine. 
C \ Results. The predictions of our models show an encouraging agreement with the demand of the self-enrichment scenario for what concerns 
i the abundances of oxygen, aluminum, fluorine and magnesium. The question of sodium is more tricky, due to the large uncertainties of the 

cross-sections of the Ne-Na cycle. The present results show that only a relatively small range of initial masses (M = 5, 6M G ) can be responsible 

for the self enrichment. 

<N ■ 
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1 . Introduction ture of all the GCs investigated is that the C+N+O sum is con- 
stant within a factor ~ 2 (Ivans et al. 1999; Cohen & Melendez 

^ , Deep spectroscopic investigations in the last decades have 2005) 

shown that Globular Clusters (GC) stars are not chemically ^ , .. u , .• . tl 

v ' 3 Two explanations, possibly acting together, are currently 

homogeneous samples, rather they show clear trends involving • .. . , . , • , , • , ,. N 

fe F ' 3 6 investigated to explain such chemical anomalies, a) the 

the chemical abundances of some light elements (Kraft 1994), , • • • /T -,r>o\ i • i 

& v " deep mixing scenario (DPS): a large non canonical extra- 

like magnesium, aluminum, oxygen, fluorine and sodium; the • • . , nr ,r, • u » u ■ j 

fe ' > je > mixing, while the stars ascend the RGB, might push inwards 

same behaviour is not followed by halo field stars, which .. , . , , , , 

3 the convective envelope to reach layers where a very ad- 
present star to star variations limited to carbon and nitrogen, , , • , t ,, ^,.,0 , • .. • , . , . ■ 

F fe vanced nucleosynthesis (full CNO burning) might have acted 

i.e. those elements whose surface abundances are expected to /r . ■ , „ . ^ nrv1 . , „ , , , 

r (Demssenkov & Weiss 2001; Demssenkov & Vandenbergh 

change following the first dredge-up and the mixing possibly , mi „ . , . , , mo , , . 1£ . , 

& & 2003; Demssenkov et al. 1998); b) the self-enrichment sce- 

following the bump on the red giant branch (hereinafter RGB). . , C r-o\ 1 c \r^-n ^ ■ 

b r b v ' nario (SES): an early generation of massive AGBs, evolving 

In almost the totality of the GCs investigated, is present a main • • . . . . . . „ ■• . . , , • 

3 in the GC, eject into the interstellar medium material which is 
stellar population, whose surface chemistry is in agreement , , , .. , .■ f .■ f , 

F F 3 & nuclearly processed, stimulating the formation of a second gen- 

with the standard a-enhanced abundances, and a second com- .. f . , , • . , ■ . 

eration of stars, whose chemistry would trace the composition 
ponent, whose surface abundances of the afore mentioned ele- f e ,. ,,, e , „ „ t-. ^ . 

F of the medium from which they formed (Cottrell & Da Costa 

ments define well determined patterns (for a recent update, see 1931 ■ Ventura et al 2001)0. 

Carretta 2006; Smith et al. 2005). These stars show depletion rr,', c tU m[ , . c „ , . , , 

' r The capability of the DPS to fully explain the observed 

of oxygen and fluorine, while the abundances of sodium and , ■ • , ■ • ■ , ■ . . • 

3 fe abundance patterns was seriously undermined by the detection 

aluminum are enhanced with respect to the solar values; the 

question of a possible depletion of magnesium by ~ 0.2 dex is 1 Here we do not consider the alternative SES proposed by Maeder 

still under debate (Cohen & Melendez 2005). A common fea- & Meynet (2006), Prantzos & Charbonnel (2006) and described by 

Decressin et al. (2007), namely that self-enrichment is produced by 
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of the same chemical anomalies even in main sequence and 
sub-giant stars of some GCs (Gratton et al. 2001), because the 
interior of these stars is not expected to reach temperatures so 
high to ignite a very advanced nucleosynthesis. 

The question concerning the SES is still open, basically 
because many ingredients used in the stellar evolution theory, 
whose physical formulation is not directly based on first prin- 
ciples (e.g. convection, mass loss rate, extra-mixing), have a 
strong impact on the physical and chemical aspects of the AGB 
evolution (Ventura & D'Antona 2005a;b); also, for one of the 
elements involved in the observed trends, i.e. sodium, the un- 
certainties associated to the relevant cross-sections (~ 3 dex, 
Hale et al. 2002; 2004) are such to render the results obtained 
very uncertain (Ventura & D'Antona 2006). Within the con- 
text of the Mixing Length Theory (MLT) modeling of turbu- 
lent convection (Vitense 1953), the theoretical ejecta of AGBs 
point against the SES hypothesis, because the great number 
of third dredge up (TDU) episodes, associated with a modest 
nucleosynthesis achieved at the bottom of the outer convec- 
tive envelope, favour a large growth of the overall C+N+O 
abundances, and inhibits the depletion of oxygen (Fenner et 
al. 2004; Denissenkov & Herwig 2003). On the other hand, 
Ventura & D'Antona (2005b), studying the effects of chang- 
ing the treatment of convection on the AGB modeling, show 
that when the Full Spectrum of Turbulence (FST, Canuto & 
Mazzitelli 1991) is used to model the convective regions, the 
FST high efficiency of convection favours higher tempera- 
tures at the bottom of the convective envelope, thus an ad- 
vanced nucleosynthesis (which is usually referred to as "Hot 
Bottom Burning", HBB), a larger luminosity, and consequently 
a smaller number of TDU episodes, thus keeping the C+N+O 
sum almost constant. 

An appealing prediction of the SES, and more generally of 
the role which AGBs may play in the pollution of the interstel- 
lar medium within GCs, is that their yields are expected to be 
helium rich, as a consequence of the deep second dredge up 
(hereinafter II DUP) experienced particularly by the most mas- 
sive models. A helium content Y> 0.30, as found in previous 
investigations (Ventura et al. 2001), may be at least partially 
consistent with the existence of helium rich stellar populations 
in GCs, which allow the interpretation of the morphology of ex- 
tended horizontal branches in some GCs (D'Antona et al. 2002; 
D'Antona & Caloi 2004; Lee et al. 2005); the presence of a he- 
lium rich population was also invoked to explain the presence 
of a blue main sequence in NGC 2808 (D'Antona et al. 2005b; 
Piotto et al. 2007) and wCen (Bedin et al. 2004; Piotto et al. 
2005). 

The possibility that DPS and SES might act together in or- 
der to explain the O-Na anticorrelation observed in a few giants 
of Ml 3 has recently been explored by D'Antona & Ventura 
(2007), and was previously suggested by Denissenkov et al. 
(1998) to account for the Mg-Al trend in the same cluster. 

Any SES model faces with a serious problem: it is 
difficult to understand how is it possible that the self- 
enriched population is very abundant in most of the 
examined clusters — see, e.g., |D'Antona & Caloi(2007)l 
concerning the ratio of the normal to self-enriched popu- 
lations derived from the analysis of the Horizontal Branch 



(HB) morphology. In part icular , for NGC 28 08, both the 
HB - [D'Antona & Caloi(2004)l |D'Antona et al.(2005b)r- 
the main sequence splittings — [D'Antona et al.(2005b)[ 
Piotto et al.(2007)| — and the spectroscopic evidence — 
Carretta et al.(2006)| — indicate that about half of the 
cluster stars is self-enriched. In order to accomplish 
this we need that: 1) either the initial mass function 
(IMF) of the first generation stars is highly anomalous, 
and is peaked at the intermediate mass stars; 2) or the 
IMF is more or less normal, the initial cluster mass is 
much larger, by about a factor 10, than the final mass, 
and the stars of the first generation have been preferen- 
tially lost, as discussed, e.g., in |D'Antona & Caloi(2004)| 
and |Prantzos & Charbonnel(2006)l In addition, 



|Bekki et al.(2007)| suggest that all GCs formed within 
dwarf galaxies, so that the cluster formation may take ad- 
vantage of the infall from all the gas lost by massive AGBs 
evolving in the galaxy and falling into the protocluster 
potential well. 

In this paper we model the AGB evolution by releasing 
some of the approximations of our previous works. We fo- 
cus on an intermediate metallicity ([Fe/H]~ -1.3). We post- 
pone to forthcoming papers the analysis concerning more metal 
rich and metal poor composition. Our approach is as follows: 
we present and discuss our intermediate mass models with 
metallicity Z=0.001, which were calculated with the latest and 
most updated physical inputs. Consistently with the hypothe- 
sis which we want to test here, we adopt an a-enhanced mix- 
ture. We compare the yields of the individual elements with the 
abundances observed in the GCs having metallicity appropriate 
for the present investigation, taking care to restrict our compar- 
ison to those stars which are scarcely evolved, to rule out any 
possible contamination from canonical or non-canonical extra- 
mixing which the stars may experience during the RGB evo- 
lution. The effects of the uncertainties related to mass loss and 
the nuclear reaction rates of the Ne-Na cycle and the Mg-Al 
chain are also examined. 



2. The ATON stellar evolution code 

The stellar evolution code used in this work is ATON3. 1 ; a full 
description of the numerical structure can be found in Ventura 
et al. (1998). Here we briefly recall the main physics inputs, 
with the most recent updates. 

2.1. Opacities 

We adopt the latest opacities by Alexander & Ferguson (2005) 
at temperatures smaller than 10000 K and the OPAL opaci- 
ties in the version documented by Iglesias & Rogers (1996). 
For both the above treatments we have the choice between the 
opacities corresponding to solar mixture (Grevesse & Sauval 
1998), and to an a-enhanced mixture [a/Fe] = 0.4. The con- 
ductive opacities are taken from Poteckhin (2006, see the web 
page www.ioffe.rssi.ru/astro/conduct/), and are harmonically 
added to the radiative opacities. 
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2.2. Equation of state 

Tables of the equation of state are generated in the (gas) 
pressure-temperature plane, according to the method described 
in D'Antona et al. (2005a). For mixtures including hydrogen, 
we use the latest OPAL EOS (2005), overwritten in the pres- 
sure ionization regime by the EOS by Saumon, Chabrier & Van 
Horn (1995). The EOS is extended to the high density, high 
temperature domain according to the treatment by Stoltzmann 
& Blocker (2000). After these large tables are written, for a 
given Z , six values of each physical quantity are computed for 
six different X values. A cubic unidimensional spline provides 
the interpolation for any input value of X. The six tables for 
H/He and given Z are supplemented by 15 tables of He/C/O in 
which the EOS is directly computed according to Stolzmann 
& Blocker (2000) as the non-ionized regions are not present 
in stellar structure following helium ignition. The interpolation 
among the 15 tables is performed using triangles in the plane 
C/O, as the stechiometric condition is Y = 1 — Xc - Xo- 

2.3. Convection 

The thermodynamic description of the regions unstable to con- 
vective motions can be addressed either within the context of 
the traditional MET formulation (Vitense 1953), or by the FST 
model (Canuto et al. 1996). 

2.4. Mixing 

Mixing of chemicals within convective zones can be addressed 
within the instantaneous mixing framework or by a diffusive 
approach. In this case, for each chemical species a diffusive- 
like equation (Cloutman & Eoll 1976) is solved: 

dX; i dX, \ d r i *> dXi n 

dt v ot 'nvd dm r l dm,. 1 

where D is the diffusion coefficient, for which, given the con- 
vective velocity v and the scale of mixing I, a local approxima- 
tion (D ~ |vf) is adopted. 

The borders of the convective regions are fixed according 
to the Schwarzschild criterium. It is also possible to consider 
extra mixing, by allowing convective velocities to decay ex- 
ponentially from the formal border, with an e-folding distance 
described by the free-parameter ( (see Ventura et al. (1998) for 
a complete discussion regarding the variation of convective ve- 
locities in the proximities of the convective borders). 

2.5. Mass loss 

Mass loss can be treated according to different prescriptions. 
It is possible to adopt the classic Reimers' treatment, the 
Vassiliadis & Wood (1993) formulation, or the prescription 
given by Blocker (1995). In this latter case the strong increase 
of the mass loss rate during the AGB evolution is modeled by 
multiplying the Reimers' rate by an ad hoc luminosity power: 
the final expression is 

M = 4.83 x lQ- 22 r] R M- 3 l L 3J R (2) 



where tj r is the free parameter entering the Reimers' prescrip- 
tion. 

2.6. Nuclear network 

The nuclear network includes 30 elements (up to 31 P) and 64 
reactions. The full list of the 30 chemicals and of the reactions 
included can be found in Ventura & D'Antona (2005a). 

The relevant cross section are taken from the NACRE com- 
pilation (Angulo et al. 1999), with only the following excep- 
tions: 

1. 14 N(p,y) 15 (Formicola et al. 2004) 

2. 22 Ne(p,y) 23 Na (Hale et al. 2002) 

3. 23 Na(p,y) 24 Mg (Hale et al. 2004) 

4. 23 Na(p,a) 20 Ne (Hale et al. 2004) 

3. The results of the most recent spectroscopic 
surveys of intermediate metallicity GC 

Several clusters of intermediate metallicity (i.e. -1.5 < 
[Fe/H]< -1.2) have been investigated spectroscopically in the 
last decade: among these, the most extensively studied are 
M3,M4,M5,M13,NGC6752 and NGC6218. The first analyses 
were focused on bright giants, and only in past few years turn- 
off (TO) and sub-giant branch (SGB) stars were studied. Bright 
giants are of little help for the main target of this paper, since 
the abundances observed might be the result of some mix- 
ing episode from the bottom of the convective envelope tak- 
ing place during the RGB evolution; on the other hand, TO 
and SGB stars never reach internal temperatures sufficiently 
high to trigger any advanced nuclesynthesis, so any chemical 
anomaly must have been imprinted in the matter from which 
they formed. Studies focused on the TO and SGB stars are 
therefore the most useful for the present investigation. 

In the analysis of the observational results, we will consider 
the uncertainties associated to the data, and the constant offsets 
among the different observers, mainly due to differences in the 
assumed solar abundances and in the temperature scale. 

The oxygen-sodium and magnesium-aluminum anticorre- 
lations are by far the most studied trends. For NGC 6752, 
the works related to TO and SGB stars (Gratton et al. 2001; 
Carretta et al. 2005) and to bump stars (Gratton et al. 2005), 
evidentiated a well defined stellar population in which the oxy- 
gen was depleted up to 0.8 dex ([0/Fe]~ -0.4) compared to 
the "standard" abundance (fO/Fe]~ +0.4), while sodium was 
enhanced at most by 0.4 dex. In Gratton et al. (2001) it is 
also shown that aluminum might be antic orrelated with magne- 
sium, the most peculiar stars showing an aluminum enhance- 
ment slightly lower than 1 dex, and a magnesium spread of 
~ 0.3dex. For this cluster Pasquini et al. (2005) studied the 
lithium content of TO stars, showing a clear anticorrelation 
between sodium and lithium, the sodium rich stars having a 
lithium content A(Li)=^l a factor of 2 smaller than the value 
detected in the stars with standard sodium. 

The O-Na anticorrelation was recently confirmed in the 
investigation of NGC 6218 performed by the same group 

2 We use the notation A(Li)=log(Li/H)+12 (in number) 
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(Carretta et al. 2006): for the stars below the RGB bump, the 
authors find a maximum extent of the oxygen depletion by 
~ ldex, and a maximum sodium enhancement of 0.5 dex. 

Ramirez & Cohen (2003) analyzed spectroscopically stars 
belonging to M5. For the stars near the TO, they found a stel- 
lar population where sodium was enhanced up to [Na/Fe]=0.3 
dex, while for oxygen they could only find upper limits for the 
estimated abundances: this is the main reason why the authors 
could not detect any clear Na-O anticorrelation. 

The case of M4 is hard to analyse, as only red giants spectra 
are available. A first extensive exploration by Ivans et al. (1999) 
evidentiated an anticorrelation among the CN weak stars, with 
a depletion of oxygen by 6 [0/Fe]~ 0.6dex, and an enhance- 
ment in sodium up to [Na/Fe]=0.3 dex. The most oxygen poor 
stars in this subsample are aluminum rich ([Al/Fe]=0.8 dex). 
The magnesium abundance is constant within 0.2 dex. A more 
recent work by Smith et al. (2005), unfortunately restricted to 
giants, confirmed the extent of the oxygen depletion (keeping 
into account a clear constant offset, which may be understood 
in the left-upper panel of their fig. 4) and of the aluminum en- 
hancement found by Ivans et al. (1999); they also found a clear 
anticorrelation between sodium and fluorine, the most sodium 
rich stars being depleted in fluorine by [F/Fe]=-0.8 dex. 

M3 and Ml 3 have been extensively studied in the literature. 
Sneden et al. (2004) evidentiated a clear O-Na anticorrelation, 
which, for the high gravity stars (i.e. those well below the RGB 

bump) extends to [O/Fe] 0.4, and [Na/Fe]=+0.4. They also 

confirm, even for this high gravity population, that the most 
oxygen poor and sodium rich stars are largely aluminum en- 
hanced, with [Al/Fe] = l. The authors also claim the detection 
in both clusters of a magnesium - aluminum anticorrelation, 
but this result has been recently argued by Cohen & Melendez 
(2005), who confirm the extent of the oxygen depletion found 
by Sneden et al. (2004), but also limit the sodium enhancement 
for the stars below the RGB bump to [Na/Fe]=+0.3. 

On the basis of the results quoted above, we may sum- 
marize the following chemical features which characterize the 
stars in the GCs which should have not experienced any "in 
situ" mechanism to change their surface abundances: 

1. The oxygen depletion is at most by 5[0/Fe] ~ 0.8 

2. The most aluminum rich stars show up an 27 A1 abundance 
of[Al/Fe]=l 

3. Sodium is anticorrelated to oxygen. The extent of the 
largest sodium enhancement is [Na/Fe]=0.4, but this result 
is uncertain by at least 0. 1 dex 

Also, we may mention other two points, which must still be 
confirmed by detailed spectroscopic investigations: 

1 . Magnesium might be correlated to oxygen, and anticorre- 
lated to sodium and aluminum 

2. Fluorine might be strongly anticorrelated to sodium, but 
this result was found at the moment only for giant stars 

3. Lithium is correlated to oxygen: the oxygen poor popula- 
tion shows a lithium content smaller by a factor 2 



4. The physical properties of the intermediate 
mass models 

The models were evolved from the pre-MS through the whole 
AGB phase. The chemistry adopted is typical of the interme- 
diate metallicity GCs, i.e. Z=0.001 and Y=0.24. The mixture, 
consistently with the main target of this investigation, is as- 
sumed to be ff-enhanced, with [a/Fe]=+0.4; the individual 
abundances are taken from Grevesse & Sauval (1998). The 
range of masses involved is 3M© < M < 6.3M Q ; the limits 
were chosen to restrict the analysis to stars achieving HBB 
during the AGB phase, and not undergoing any carbon igni- 
tion in the interior. Convection was modelled according to the 
FST treatment. Due to the importance that CNO burning within 
the most internal part of the convective envelope may have on 
the AGB evolution (Herwig 2005; Mazzitelli et al 1999), we 
adopted the diffusive treatment in all the models presented here. 
In the phases preceeding the AGB evolution, a free parameter 
( = 0.02 was used to model the exponential decay of velocities 
within regions stable against convection; this is in agreement 
with the calibration in Ventura et al. (1998). For the whole AGB 
phase no extra-mixing was assumed from any convective bor- 
der. Mass loss was modelled according to Blocker (1995), with 
the free parameter entering eq.2 set to t/r = 0.02, according to 
the calibration given in Ventura et al. (2000). 

The main physical properties of the models are reported in 
Table Q] Cols. 2 to 4 show the duration of the core H- and He- 
burning phases (Myr) and of the AGB phase (Kyr). 

Within stars of intermediate mass, after the extinction of 
the CNO burning shell following the core He-burning phase, 
the bottom of the convective envelope sinks inwards, down to 
layers precedently touched by nuclear burning; this is the sec- 
ond dredge-up (Iben 1991) Col. 5 reports the total mass previ- 
ously involved in CNO burning mixed with the convective en- 
velope during the II DUP; in the most massive models, where 
the amount of dredged-up matter is higher, the surface chem- 
istry is altered, with a decrease of the surface oxygen abun- 
dance and an increase of the sodium and helium mass fractions. 

These results are consistent with the old models calculated 
with the same metallicity, presented in Ventura & D'Antona 
(2005b, see Table 1). In cols. 6-10 of Table Q] we give some 
details of the AGB evolution, namely the core mass at the be- 
ginning of the Thermal Pulses (TPs) phase, the maximum lu- 
minosity and the maximum temperature reached at the bottom 
of the convective zone (T\, ce ), the number of the first TP fol- 
lowed by a TDU episode, and the maximum efficiency of the 
TDU during the whole stellar AGB life0. 

FigQ]shows the variation with mass (decreasing during the 
evolution) of the luminosity and 7b ce in some of our models. 
We note in all cases, with the only possible exception of the 
3M model, a rapid increase of the luminosity after the very 
first TPs, associated to an increase of T\, ce . This is a conse- 
quence of the high efficiency of the convective model adopted, 
as correctly predicted by Blocker & Schonberner (1991), and 

3 We use here the usual descriprion of the efficiency of the TDU 
in the terms of the quantity A, defined as the ratio of the total mat- 
ter dredge-up after a TP to the mass by which the H-exhausted core 
advanced from the previous TP 
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Fig. 1. Variation with the total mass of the luminosity (Left) and Temperature at the bottom of the convective envelope (Right) 
of some intermediate mass models during the AGB phase. For clarity reasons, only models corresponding to initial masses 
3,4,5,6M are shown. 



later confirmed by D'Antona & Mazzitelli (1996) and Ventura 
et al. (2000). A detailed comparison of the results obtained with 
various efficiencies of the convective model can be found in 
Ventura & D'Antona (2005b). 

The large luminosities attained by our models have two im- 
portant consequences: 

1. A fast decrease of the mass of the envelope, with a conse- 
quent small number of TPs, and therefore of TDUs 

2. A very advanced nucleosynthesis at the bottom of the en- 
velope 

In the most massive models, mass loss is so large that they 
reach their maximum luminosity after a few TPs; they loose 
their envelope so rapidly, that TDU takes place only in the 
latest stages of their AGB evolution, and with a very modest 
efficiency (see col. 10 of Table Q}. We stress that the larger is 
the luminosity, the faster is the general cooling of the structure 
determined by the gradual loss of the mass of the envelope: 
the maximum temperature achieved at the bottom of the ex- 
ternal convective zone reaches a maximum asymptotic value 
(see col. 8 of Table [T), which, for large M, turns out to be in- 
dependent of the total initial mass of the star; in the present 
computations, this upper limit is T'bce = 110MK. 

5. The chemistry of the ejecta 

For each model, we calculate the average mass fractions in the 
ejecta, for the chemical elements included in our network. The 
results for the species of interest for this work, are presented 
in Table [2] For any isotope A, we give the quantity [A/Fe], de- 
fined as [A/Fe] = log(A/Fe) e jecta - log(A/Fe) ; the abundances 
are mass fractions. The only exceptions are helium and lithium, 
for which we list, respectively, the mass fraction Y and the stan- 
dard A(Li) quantity. Note that Mg in col. 9 refers to the total 
magnesium abundance, and R(CNO) in the last column gives 
the ratio between the total C+N+O in the ejecta and the initial 
value. 

5.1. The oxygen depletion 

The depletion of the surface oxygen abundance requires a 
strong HBB at the bottom of the convective envelope, as the ac- 



tivation of the full CNO cycle demands temperatures approach- 
ing 100MK, which, within the context of the AGB modelling, 
are attainable only via a very efficient description of convec- 
tion (Ventura & D'Antona 2005b). The present models have 
been calculated with the FST prescription, thus, at least in the 
most massive models, we expect to reach such a high T'bce, as 
can be seen in the right panel of Fig[T] 

The left panel of Figj2] shows the variation with the total 
mass of the surface oxygen abundance for the masses examined 
here, with the only exception of the 6.3M model, which was 
omitted for clarity reasons, being very similar to the 6M case. 
We see that oxygen burning starts efficiently shortly after the 
beginning of the AGB phase in all the models more massive 
than 4M . In the same panel we may easily note the effects of 
the TDU, which, when sufficiently efficient, increases the oxy- 
gen content of the envelope, because carbon and oxygen rich 
material is dredged-up from the ashes of the precedent 3a burn- 
ing shell which forms during each TP. A strong oxygen deple- 
tion is thus inhibited by repeated and efficient TDU episodes, 
so that in the less massive models, which undergo many TDUs, 
and hardly reach the temperatures requested to ignite oxygen 
burning, the oxygen is indeed produced rather than destroyed 
within their envelopes. With increasing mass, we shift progres- 
sively to a situation where the final oxygen abundance is a del- 
icate compromise between the depletion triggered by the acti- 
vation of the full CNO cycle and the increase due to the TDU 
(4.5M < M < 5.5M ), to end up with the most massive models 
(M >6M Q ), in which oxygen can be eventually depleted by a 
factor ~ 20 compared to the initial value (see the track corre- 
sponding to the 6M model in the left panel of FigO. 

In the right panel of Figj2] we show the oxygen content of 
the ejecta of the models calculated, in terms of [O/Fe], to al- 
low a more straight comparison with the observations outlined 
in Sect. 2. We note a very high oxygen content in the ejecta 
of the models with masses M < 4M , consistently with the 
previous discussion; the model with initial mass 4M shows an 
oxygen content unchanged compared to the initial a-enhanced 
value. The mass-oxygen trend is progressively decreasing with 

mass, and reaches a plateau value of [O/Fe] 0.4dex for all 

the masses M > 5M . The reasons why this lower limit ex- 
ists is twofold: a) on the one hand (see col. 8 in Table [U, we 
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Table 2. Chemical composition of the ejecta of intermediate mass models 



M/M Q 




Y 


A(Li) 


[ 12 C/Fe] 


[ 14 N/Fe] 


[ 16 0/Fe] 


[ 19 F/Fe] 


[ 2J Na/Fe] 


[Mg/Fe] 


[ 27 A1/Fe] 


R(CNO) 


3.0 


0.02 


.248 


2.77 


0.84 


2.21 


0.92 


0.10 


1.16 


0.57 


0.65 


9.6 


3.5 


0.02 


.265 


2.43 


0.51 


2.18 


0.77 


-0.26 


1.30 


0.55 


0.66 


7.9 


4.0 


0.02 


.281 


2.20 


0.14 


2.02 


0.44 


-0.61 


1.18 


0.48 


0.55 


4.9 


4.5 


0.02 


.310 


2.00 


0.12 


1.89 


0.19 


-0.90 


0.97 


0.43 


0.85 


3.1 


5.0 


0.02 


.324 


1.98 


0.13 


1.70 


-0.06 


-1.16 


0.60 


0.35 


1.02 


2.1 


5.5 


0.02 


.334 


1.93 


-0.41 


1.51 


-0.35 


-1.39 


0.37 


0.28 


1.10 


1.3 


6.0 


0.02 


.343 


2.02 


-0.62 


1.35 


-0.40 


-1.36 


0.31 


0.29 


1.04 


0.97 


6.3 


0.02 


.348 


2.06 


-0.68 


1.33 


-0.37 


-1.28 


0.30 


0.30 


0.99 


0.94 


5.0 


0.01 


.327 


1.79 


0.00 


1.83 


-0.14 


-1.30 


0.67 


0.40 


1.20 


2.77 


5.0 


0.04 


.323 


2.49 


-0.37 


1.58 


-0.05 


-1.00 


0.70 


0.39 


0.80 


1.71 


6.0 


0.01 


.345 


1.83 


-0.40 


1.49 


-0.42 


-1.45 


0.27 


0.22 


1.07 


1.29 



Fig. 2. Left: Variation with the total mass of the surface oxygen mass fraction of the AGB models with masses 3M Q < M < 6M ; 
the strong depletion in the most massive models is a clear signature of strong HBB at the bottom of the convective envelope. 
Right: oxygen content of the ejecta as a function of the initial mass; a plateau value of [O/Fe]=-0.4 is reached for the highest 

masses 



saw that there is an upper limit for T\, ce , which therefore limits 
the degree of O-burning which may be achieved at the bottom 
of the envelope; b) the luminosity in the most massive mod- 
els is so high that they loose mass rapidly, already from the 
very first TPs, when the oxygen abundance is still large (see 
the different slope of the Oxygen-Mass relation characterizing 
the 6M Q model compared to the other masses in the left panel 
of Figf2|. We will show that changing the mass loss descrip- 
tion does not change substantially this conclusion. Finally, we 
note that this lower limit for [0/Fe]is in good agreement with 
the lowest oxygen abundances measured in TO and SGB stars 
belonging to intermediate metallicity GCs. 

5.2. Aluminum production and the activation of the 
Mg-AI cycle 

Cols. 10 and 11 of Table [2] report the Mg and Al content of 
the ejecta. While the observed spread in the Mg abundance, as 
outlined in Sect.3, is still debated, and in any case restricted 
to 0.2-0.4 dex, the stars in GCs showing up the largest degree 
of oxygen depletion are also strongly enriched in aluminum, 
with [Al/Fe] = l in the most extreme cases (Gratton et al.2001; 
Sneden et al. 2004). 

In the AGB modelling, an aluminum production at the 
bottom of the envelope is made possible by the activation 
of the Mg-AI chain, which, favouring proton captures by the 
heavy isotopes of magnesium, eventually leads to the synthe- 
sis of 27 Al (Denissenkov & Herwig 2003; Ventura & D' Antona 
2005a). Also TDU plays a role which is not negligible, as a 
deep penetration of the convective envelope within the region 
precedently touched by helium burning may bring to the sur- 
face 25 Mg and 26 Mg synthetized via a- captures by 22 Ne nu- 
clei; these isotopes, once the CNO burning shell is reactivated, 
undergo proton capture and produce aluminum. 

Figj3] shows the evolution of the surface aluminum content 
in the models presented here. We note in all cases a trend in- 
creasing during the evolution; in the most massive models, due 



Fig. 3. The variation during the AGB phase of the surface alu- 
minum abundance, for the same models reported in figH] Note 
the combined effects of HBB and of TDU to increase the sur- 
face aluminum content 

to the stronger nucleosynthesis activated, the aluminum pro- 
duction is larger. Even in this case, as it was for the depletion of 
oxygen, we note an upper limit to the aluminum enhancement, 
which can also be seen in col. 10 of Table|2]to be [Al/Fe] = l. 
The reasons for the impossibility of a larger Al production are 
the upper limit of T\, ce = 1 10MK found in our models, the neg- 
ligible extent of TDU in the most massive stars, preventing the 
transport of the magnesium isotopes from the interior to the 
surface, and the strong mass loss suffered from these models 
at the beginning of the AGB phase, which favours the ejection 
into the ISM of material which is not aluminum rich. 

We see from col. 1 1 of Table[2]that the ejecta of all the mod- 
els are characterized by 0.5 <[A1/Fe]< 1, and are therefore 
fully consistent with the aluminum content of the stars with the 
anomalous chemistry outlined by the spectrospic investigations 
presented in Sect.3. 

The total magnesium abundance, as it is evident from the 
previous discussion, is the result of the balance between the in- 
crease of Mg determined by the TDU, and the depletion due to 
proton captures during HBB. This explains the negative trend 
with mass which can be seen in col. 9 of TablefJ] 

We underline here the striking difference between our find- 
ings and the results obtained by Fenner et al. (2004, see the 
bottom panel of their fig. 1), where they found that the most Al- 
rich stars were also magnesium rich. This was a result of the 
effects of many TDUs, enriching the envelope with the heavy 
magnesium isotopes produced in the 3a shell; it is the different 
treatment of convection between the two sets of models leading 
to this discrepancy, because the use of the FST model reduces 
the number of TPs and TDUs. 

Contrary to oxygen, the predictive power of our results for 
Al is undermined by the range of uncertainties related to the 
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Fig. 4. Variation with mass of the surface sodium abundance 
within our standard models 

Fig. 5. Variation of the sodium surface content within models 
with initial masses 4,5,6M calculated with the recommended 
values of the cross-sections of the Ne-Na cycle (dotted), and 
those maximising sodium production (see text for details). The 
dashed tracks show the results obtained when the upper limits 
of the a-capture reactions by 22 Ne nuclei are adopted. 

relevant cross-sections. Izzard et al. (2007) evidentiated that in 
massive AGBs the yields of 27 Al is affected by the uncertainties 
connected to both the 26 Mg(p,y) 27 Al and 26 Al(p,y) 27 Si reaction 
rates: we used the upper NACRE limit for these reactions in 
the present investigation, but we keep in mind that these results 
have an associated uncertainty which may be estimated to be 
around 0.3-0.4 dex. 

5.3. The puzzling behaviour of sodium 

The debate regarding the amount of sodium which may be syn- 
thesized within AGBs is still open, due to the large uncertain- 
ties associated to a) the cross-sections of the reactions involved 
in the Ne-Na cycle; b) the cross sections of the <r-captures by 
22 Ne; c) other physical inputs which play a role in determining 
the sodium content within the envelope of these stars. 

As outlined by Ventura & D'Antona (2006), the surface 
sodium abundance first increases due to the II DUP, then, par- 
ticularly in AGB models calculated with an efficient convec- 
tive model, it is further produced by burning of the dredged-up 
22 Ne, and is later decreased when the rate of destruction ex- 
ceeds that of production. Any TDU favours sodium production, 
due to dredging up of primary 22 Ne synthesized via a-capture 
within the convective shell which forms during the TP. This 
behaviour is confirmed by FigJU which shows the variation of 
the surface sodium abundance of the evolutionary models. The 
behaviour of the M< 4M Q stars is in qualitative agreement with 
the AGB models used by Fenner et al.(2004) (see the upper 
panel of their Fig.l): we note a great increase of the sodium 
abundance, due to the dredge-up of 22 Ne which is later con- 
verted to sodium. Contrary to their findings, our more massive 
models show an opposite behaviour, because the FST convec- 
tive model favours a more advanced nucleosynthesis, with a 
partial destruction of the sodium precedently created; also, we 
recall that the higher mass loss favours a smaller number of 
TPs, thus acting against sodium production. 

The precedent discussion explains the clearly negative 
trend with mass of the [Na/Fe] values in the 9th column of 
Table [2] We focus our attention on the most massive models, 
which produce yields which we saw to be aluminum rich and 
oxygen poor: the sodium content of their ejecta is in the range 
0.3-0.4 dex, which is consistent with the sodium abundances 
derived by most of the research groups for the stars with the 
most anomalous chemical composition. 

The uncertainties connected to the cross sections have 
a dramatic impact on the value of [Na/Fe] of the ejecta, 



even more than we saw for aluminum. Concerning HBB, 
the main problems are associated to the cross section of the 
22 Ne(p,y) 23 Na reaction, which is uncertain by a factor ~ 2000 
(Hale et al. 2002); even the other reaction relevant to determine 
the correct sodium equilibrium value, i.e. 23 Na(p,a) 20 Ne, has a 
margin of uncertainty, which is however smaller (~ 30%; Hale 
et al. 2004). 

We ran three models with initial masses 4,5, 6M , where 
we used the lower limit for the 22 Ne(p,y) 23 Na reaction, and 
the upper limit for 23 Na(p,a) 20 Ne. The comparison between 
these simulations (dotted tracks) and those described in Table 
[2] (solid) is shown in FigO We note that when the cross sec- 
tions minimizing the sodium production are used, with the 
only exception of the II DUP, sodium is destroyed when T\, ce 
becomes sufficiently large to ignite proton capture by 23 Na 
nuclei. The sodium contents which we get are considerably 
lower, i.e. [Na/Fe] =0.6 for the 4M model, [Na/Fe] =0.0 in the 
M=5M case, and [Na/Fe] =-0.2 for M=6M ; these values are 
smaller than those reported in Table|2]by 0.6 dex. Note that this 
is not proportional to the reduction factor of the 22 Ne(p,y) 23 Na 
reaction, because, as pointed out by Izzard et al. (2007), once 
22 Ne is destroyed at the bottom of the envelope, no further 
sodium can be created, despite the use of a cross section for 
22 Ne burning which is a factor ~ 2000 higher. This is also con- 
sistent with the investigation by Ventura & D'Antona (2006), 
who pointed out that a decrease of the 23 Na(p,a) 20 Ne reaction 
by a factor 2 could reconcile better the sodium content of the 
ejecta of the most massive AGBs of intermediate metallicity 
with the spectrospic measurements of GC stars. 

For the models experiencing TDU, another source of un- 
certainty for the sodium yield is provided by the cross sections 
of the a- capture reactions by 22 Ne nuclei, whose upper limit 
is 3 orders of magnitude higher than the recommended values 
(Angulo et al. 1999). These reactions determine the 22 Ne con- 
tent in the convective shell which forms during the TP, hence 
the amount of 22 Ne which may be dredged-up in the after- 
pulse phase, and therefore the quantity of sodium which may 
be sinthesized at the bottom of the envelope via proton capture 
by 22 Ne nuclei. 

The dashed tracks in Fig|5]show the results of the surface 
sodium abundance when these upper limits are used: for the 
4 and 5M the overall sodium abundance is clearly reduced. 
Differently from the previous case, the associated uncertain- 
ties on the sodium yields are not uniform with mass, but rather 
show a decreasing trend, ranging from a null effect for the most 
massive models, to 5[Na/Fe] ~ 0.15dex for M = 5M , up to 
5[Na/F <?] ~ 0.3dex for M = 4M ; these results are actually 
not surprising, since dredging up of 22 Ne becomes progres- 
sively more important in determining the sodium abundance 
the smaller is the stellar mass. 

The yields of sodium from these sources are still highly 
uncertain; only a more solid estimate of the relevant cross- 
sections may help increasing the reliability of these investiga- 
tions. 
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Fig. 6. Left: Variation as a function of the total stellar mass of the surface lithium content of the intermediate mass models during 
the whole AGB evolution. Right: The lithium content of the ejecta, as a function of the initial mass. 



5.4. The lithium problem 

Lithium is synthesized at the bottom of the convective envelope 
of AGBs at temperatures exceeding 40 x 10 6 K via the Cameron 
& Fowler (1971) mechanism. Sackmann & Boothroyd (1992) 
first showed that the use of a diffusive approach was mandatory 
to describe such a delicate interplay among nuclear and mix- 
ing time scales, which could eventually lead to lithium produc- 
tion. Ventura et al. (2000) used this approach to reproduce the 
lithium vs. luminosity trend observed in the Large Magellanic 
Cloud. 

The key factor to achieve lithium production is the activa- 
tion of the 3 He(a, y) 7 Be reaction at the bottom of the convec- 
tive envelope, which is possible when T\, ce > 40 x 10 6 K. 

The left panel of Fig|6]shows that all our models attain tem- 
peratures at the bottom of the convective envelope sufficiently 
high to ignit the Cameron & Fowler mechanism; the most mas- 
sive stars undergo a rapid consumpion of the whole 3 He avail- 
able in the envelope, so that the surface lithium content, after 
reaching a maximum value A(Li)~ 4, rapidly declines to ex- 
tremely low abundances. This process becomes progressively 
slower as the mass decreases; we end up with the 3M Q model, 
which is still lithium rich at the end of its AGB evolution. This 
discusion explains the relation between the initial mass of the 
star and the lithium content of its ejecta, shown in the right 
panel of Fig|6j all the massive models have A(Li) ~ 2, in ex- 
cellent agreement with the lithium abundance of the oxygen 
poor TO stars in NGC 6752 (Pasquini et al. 2005). 

5.5. The helium enrichment 

The AGB models discussed here experience a small number of 
TPs, thus most of the helium enrichment of the envelope takes 
place during the II DUP. Since this latter is deeper the higher 
is the initial mass of the star (see col.5 of Table [U we expect 
the helium enrichment to increase with mass, as it is confirmed 
by the results reported in col. 3 of Table [2] We find a small in- 
crease (compared to the standard Big Bang value, Y=0.24) in 
the models with mass M < 4M©; the maximum enrichment, 
for the masses close to the limit for carbon ignition, is Y=0.35. 
We discuss in a forthcoming paper (Pumo et al. 2007) the over- 
all helium enrichment due to AGB and super- AGB stars of in- 
termediate metallicity, and the robustness of these predictions, 
related essentially to the efficiency of the II DUP. 

5.6. Comparison with previous models 

Ventura et al. (2002) presented models for the evolution of stars 
of intermediate mass at various metallicities, ranging from Z = 
2 x 1 0~ 4 to Z = 0.0 1 . We compare the yields of the present work 
with the results of that investigation, for the metallicity Z = 
0.001. First, we note that the helium content is sistematically 
higher here by ~ 0.04, as can be seen by comparing the third 
column of Table [2] with the dashed curve giving the helium- 



mass trend in fig. 4 in Ventura et al. (2002). This results may be 
understood in terms of the extra-mixing from the bottom of the 
convective envelope which was assumed in the present models, 
and which was neglected in Ventura et al. (2002). 

The combined effects of the overshooting from the enve- 
lope and the different mixture adopted (solar in Ventura et al. 
(2002), a- enhanced here), favours smaller core masses in the 
present models, so that a given model here can be compared 
with a model 0.5M o less massive in Ventura et al. (2002). 

Even with this assumption, we note differences in the CNO 
yields, the overall C+N+O contribution being higher here com- 
pared to the corresponding values in Ventura et al. (2002) (see 
fig. 6 in Ventura et al. 2002, to be compared with columns 5 
to 7 in Table |2). The reason for such discrepancy is the differ- 
ent efficiency of the TDU found in the two sets of models: we 
find a maximum efficiency of A — 0.7 here (which decreases 
to A = 0.5 in the more massive models), whereas in Ventura et 
al. (2002) A could hardly reach 0.5. This change, favoured by 
the different mixture adopted, determines the differences in the 
yields obtained. 

The detailed comparison of the other yields is rendered hard 
by the differences in the nuclear cross-sections adopted for the 
relevant reactions; Ventura et al. (2002) used the NACRE com- 
pilation, whereas here we use the most updated releases present 
in the literature. 

6. The uncertainties related to mass loss 

Mass loss plays a fundamental role in the context of AGB evo- 
lution: it is the reduction of the convective envelope via stel- 
lar winds which eventually halts the TP phase, and leads to 
PN ejection. Also, the description of mass loss determines the 
number of TPs experienced by the star, and, in turn, the number 
of TDUs (Schonberner 1979). 

Ventura et al.(2000) calibrated the parameter j]r entering 
eq.(2) by reproducing the luminosity function of lithium rich 
stars in the Magellanic Clouds; the chemistry of the stars ex- 
amined in the present work is different, thus leaving some room 
for a possible variation of It is therefore essential to under- 
stand to which extent the yields which we obtain depend on 
this choice. 

On the physical side, a larger mass loss leads the star 
evolve at lower luminosities; the degree of the nucleosynthesis 
achieved at the bottom of the convective envelope is reduced, 
because we have smaller temperatures. This is confirmed by 
figH where we show the evolution of the luminosity and of 
7bce m models with initial masses 5 and 6M , calculated with 
different 77RS. 

Chemically, the situation is complex, and not all the iso- 
topes follow the same trend with M; this is dependent on the 
modality with which any chemical species is synthesized (or 
destroyed) within the convective envelope. The last three lines 
of Table |2]report the chemistry of the ejecta of the models of 5 
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and 6M Q calculated with a different mass loss rate with respect 
to our standard case. 

The behaviour of lithium is the most linear. As can be seen 
in the left panel of fig HO lithium is produced only during the 
first TPs, so that the average lithium content of the ejecta is 
determined essentially by the mass lost by the star during this 
phase. A large M allows a larger release of lithium rich mate- 
rial during the early phases of the AGB evolution, with a conse- 
quent increase of A(Li): for a given mass, we see in Table|2]an 
almost linear trend A(Li)-?7 S . A similar, straight path is also fol- 
lowed by fluorine, which is destroyed during the very first TPs, 
so that the overall yield is determined by the strength of the 
stellar winds during the first TPs; in this case the slope of the 
[F/Fe]-?7R relation is lower compared to lithium, because the 
fluorine consumption is faster than the duration of the whole 
phase of lithium production and destruction. 

Other elements show up a less defined behaviour, because 
their yield is determined by the nucleosynthesis at the bottom 
of the envelope during the whole AGB phase, and also by the 
TDU. Interestingly, we find that in some cases the yields are 
not very sensitive to the mass loss rate adopted. A typical ex- 
ample is oxygen, for which we show the variation of the sur- 
face abundance in the right panel of figH] When mass loss is 
reduced, the tendency of the oxygen abundance of the ejecta 
to diminish (less oxygen-rich matter during the first TPs is lost 
by the star) is partly compensated by the larger number of TPs, 
which act to increase the surface oxygen content. The clear- 
est example is the 6M model, for which a smaller rju (dotted 
track) leads to a smaller surface oxygen abundance until the 
total mass of the star drops to 3M Q , but to a larger mass frac- 
tion in the latest evolutionary phases, when the TDU becomes 
efficient. As can be seen in Table 12 the oxygen content of the 
ejecta is almost independent of and so is the general con- 
clusion which we reached in Sect. 5. 1 concerning the maximum 
depletion of oxygen obtainable by these models: playing with 
mass loss leaves unchanged the maximum depletion of oxygen 
obtainable at these metallicites, leading to a minimum oxygen 
content of the ejecta [O/Fe]=-0.4. 

A behaviour similar to oxygen is also followed by sodium, 
again because its abundance in the envelope is a balance be- 
tween destruction via proton capture and production via dredg- 
ing up of 22 Ne. Independently of mass, the related uncertainty 
is (5[Na/Fe]~ 0.1 dex. 

The aluminum abundance of the ejecta decreases with M, 
because a larger mass loss favours a larger release of mass at the 
beginning of the AGB phase, when the aluminum has not yet 
been synthesized in great amounts. The smaller number of TPs 
at large M tends to increase the 27 Al abundance indirectly, via 
dredge-up of the heavy magnesium isotopes, which later form 
aluminum via proton capture. The variation of [Al/Fe]with the 
mass loss rate is [5[Al/Fe] ~ O.ldex). 

Helium is not influenced by the details of the mass loss de- 
scription, because the only substantial change of its surface 
abundance takes place during the II DUP, and remains approx- 
imately constant during the whole following AGB phase. 

Finally, we note from the last column of Table [2] that the 
CNO ratio has a negative trend with mass loss, which is a mere 



Fig. 7. Comparison between the variation with the total mass of 
the luminosity (top) and temperature at the bottom of the con- 
vective envelope (bottom) of models with initial masses 5,6M , 
calculated with different values of the parameter j]r enetering 
eq.2 

consequence of the fact that when M is high the star experience 
a smaller number of TPs. 

We may summarize the effects of the uncertainties of mass 
loss during the AGB phase on the chemistry of the ejecta as 
follows: 

1 . Oxygen is not sensitive to the mass loss rate adopted 

2. Sodium and aluminum show a positive trend with mass 
loss, but the uncertainty associated with the adopted de- 
scription of mass loss is considerably smaller than the in- 
determination due to the unknown relevant cross-sections 

3. Lithium and fluorine are more sensitive to mass loss, with 
a linear positive trend 

4. The CNO content of the ejecta is sensitive to mass loss, 
because a stronger M diminishes the number of TPs ex- 
perienced by the stars, thus favouring a smaller C+N+O 
abundance 



7. Discussion: theory vs. observations 

We discuss the self-enrichment scenario hypothesis for inter- 
mediate metallicity clusters, by asking whether the ejecta of 
the most massive AGBs can account for the chemical patterns 
traced by the abundances of the GC stars with the most anoma- 
lous chemistry. We restrict our attention on the least evolved 
stars (i.e. TO and SGB sources, or giants well below the RGB 
bump), despite the difficulties presented by their spectroscpic 
analysis, because this allows to rule out any possible change 
of the surface chemistry due to some non-canonical extra mix- 
ing while ascending the RGB; we therefore disentangle the pri- 
mordial from the evolutionary effects, and focus only on the 
abundance patterns present directly in the matter from which 
the stars formed. 

Our goal is to test the possibility that our ejecta can repro- 
duce the observed O-Na and O-Al trends, and that the stars 
showing the strongest oxygen depletion are fluorine poor and 
possibly sligthly depleted in magnesium. Finally, we compare 
our results with the recent analysis of the lithium abundances in 
NGC 6752, which indicate that oxygen poor stars deviate from 
the Spite's plateau, having a lithium content a factor 2 smaller 
than the standard value (Pasquini et al. 2005). 

In making this comparison, we keep in mind that the oxy- 
gen abundances predictions are very robust, that the aluminum 
and even more sodium mass fractions are made uncertain by 
the poor knowledge of the relevant cross-sections in the range 
of temperatures of interest here (T~ 100MK), and that lithium 
is strongly influenced by the assumed mass loss rate at the be- 
ginning of the AGB phase. 

The left panel of fig|9] shows the observed oxygen-sodium 
trend for such stars. The dashed track and the two solid lines in- 
dicate the abundances of the ejecta of our models, differing in 
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Fig. 8. Variation with the total mass of the star of the lithium (Left) and oxygen (Right) surface abundances within models of 
initial mass 5 and 6M calculated with different mass loss rates. Note the straightforward dependency of lithium on mass loss, 
compared to the more tricky behaviour of oxygen, whose destruction is first amplified by a lower mass loss rate, and later 
prevented by TDU episodes 



the chosen cross-sections for the reactions involving sodium: 
the upper solid track refers to models calculated with the up- 
per limits for the 22 Ne(p,y) 23 Na and the lower limits for the 
23 Na(p,a) 20 Ne reactions, while the lower refers to models cal- 
culated with the recommended values for the same reactions; 
the dashed track indicates the results obtained by maximizing 
the rates of the a-capture reactions by 22 Ne. 

According to our interpretation, the observed points inside 
the squared box in the right lower portion of the plane rep- 
resent stars born with the original chemistry, while those be- 
longing to the 2nd generation, whose initial chemistry traces 
the pollution by AGBs, are included within the squared box 
in the left-upper part. The remaining points, with high sodium 
and normal oxygen, may be stars formed by processed mat- 
ter mixed with remnant primordial gas (Decressin et al. 2007). 
The uncertainties related to the cross-sections strongly limit the 
predictive power of the results obtained, and the observational 
spread (£[Na/Fe]~ 0.3 for a given [O/Fe]) makes the compar- 
ison not straightforward. However, we note that the ejecta of 
models with masses M > 5M Q might account for the oxy- 
gen and sodium abundances detected in the 2nd generation of 
stars. Notice that here we are touching again the problem 
of the mass budget: if only the envelopes of stars from 5 
up to 6.3M — or at most up to 7-8M , if we can assume 
that also the superAGBs contribute with similar yields — 
can form the self-enriched stars, the gas contained in this 
mass range, for reasonable IMFs, is only a few percent of 
the total initial cluster mass. As we remarked in the in- 
troduction, we must then hypothize that the initial cluster 
mass was much larger than the present mass, and that pref- 
erentially the second generation stars have been lost during 
the long term cluster evolution. 

The right panel of Fig|9]reports the observed points in the 
O-Al plane. Even in this case we report two theoretical lines, 
obtained with different choices of the rates of the Mg-Al chain 
reactions (the upper line refers to models calculated with the 
maximum allowed values of the cross-section of the proton 
capture reactions by the two heavy magnesium isotopes). In 
this case the observed trend is well reproduced, in particular 
for the chemistry of the most oxygen poor stars, which evi- 
dentiate an aluminum enhancement by [Al/Fe]~l. In the case 
of aluminum the comparison is more straightforward, because 
the theoretical uncertainties related to the cross-sections are 
smaller. The two squares in the figure have the same meaning 
as in the left panel. 

A welcome result from this investigation is that the lithium 
content of the ejecta of the most massive models, those showing 
the strongest depletion of oxygen, is A(Li)~ 2 (see tebj2]», and 
are therefore in excellent agreement with the recent analysis of 
TO stars in NGC 6752 by Pasquini et al. (2005). Actually, the 
overall O-Li trend in the range [0/Fe]< 0.2 is well reproduced 
(see their fig. 3). 



Although the observed range of magnesium abundances is 
small, so that the existence of a real spread is still debated, our 
models indicate that a Mg-Al anticorrelation should be present, 
with a spread d>[Mg/Fe]~ 0.3dex. 

Finally, our models point in favour of a strong depletion 
of fluorine, which is burnt heavily at the very beginning of the 
AGB phase in all models whose initial mass exceeds 4M Q . This 
is in agreement with a recent investigation focused on M4 gi- 
ants (Smith et al. 2005); a robust confirm of this scenario would 
be the determination of fluorine abundances in TO and SGB 
stars. 

8. Conclusions 

We present updated model for the evolution of AGB stars of 
3M < M < 6.3M , Y=0.24 and Z=0.001. We test the self- 
enrichment scenario hypothesis for GCs of intermediate metal- 
licity, comparing the results of spectroscopic investigations of 
scarcely evovlved stars in GCs with the theoretical yields of in- 
termediate mass stars. This task demanded the computation of 
a new set of models, calculated with the latest physical updates 
concerning the equation of state and the conductive opacities, 
and with an a-enhanced mixture, to describe self-consistently 
the chemistry of the first generation of stars which form in the 
GCs. 

The high efficiency of the convective model adopted con- 
firms an important result obtained by this research group re- 
garding the AGB evolution of intermediate mass stars, i.e. the 
possibility of a strong nucleosynthesis at the bottom of the ex- 
ternal convective zone for all the masses M> 3M ; this HBB 
also favours a fast increase of the luminosity, a higher mass 
loss, and therefore reduces the number of thermal pulses ex- 
perienced by the star during the AGB phase. The present in- 
vestigation indicate that for this metallicity a maximum T\, ce = 
1 10MK may be reached for the highest masses; this upper limit, 
and the strong mass loss suffered by the most massive models, 
limits the extent of the nucleosynthesis which may be achieved 
within the envelope. 

The combination of HBB and TDU in the most massive 
models favours a strong depletion of oxygen and fluorine, a 
modest reduction of magnesium, and a large production of alu- 
minum. Sodium is also produced, via a delicate compromise 
between production by neon burning and destruction by proton 
capture. It is confirmed that lithium can be produced at the be- 
ginning of the AGB phase via the Cameron-Fowler mechanism, 
to be later destroyed due to 3 He consumption within the enve- 
lope. The matter ejected by these models is helium rich, which 
a maximum enrichment of Y=0.35 found for the 6M model. 

On the basis of the results of this work, the strongest point 
in favour of the self-enrichment scenario is the oxygen deple- 
tion, for which both the observations an the theoretical predic- 
tions indicate a maximum limit of [0/Fe] = -0.4; the theoreti- 



Ventura & D' Antona: AGB models of intermediate metallicity 



11 



Fig. 9. Left: The observed Oxygen-sodium trend in stars in GCs. The solid lines indicate the content of the ejecta of our models, 
obtained with the recommended cross-sections for the reactions involved in the Ne-Na cycle (lower track), and with the reactions 
rates maximizing sodium production (upper lines); the dashed track shows the abundances of the ejecta from the models calcu- 
lated with the enhanced cross sections for the a-captures by 22 Ne. Right: The observed Oxygen-Aluminum anticorrelation. The 
two solid lines indicate the results from our models, accroding to the choices made for the cross-section of the Mg-Al reactions. 
The observed points refer to the following works. Full triangles: NGC 6752 (Gratton et al. 2001); Open squares: M5 stars with 

V > 16 (Ramirez & Cohen 2003); Stars: M13 stars with V > 15 (Cohen & Melendez 2005); Open triangles: M3 stars with 

V > 15 (Cohen & Melendez 2005); Full points: NGC 6218 stars (Carretta et al. 2006); Open squares: High gravity M3 giants 
(Sneden et al. 2004); Full squares: High gravity M13 giants (Sneden et al. 2004) 



cal oxygen yield is robust, as it turns out to be approximately 
independent of the details of the mass loss description, and the 
relevant cross-sections are known with sufficient accuracy. 

The O-Al trend is confirmed by the present investigation, 
though the extent of the aluminum enrichment of the ejecta is 
sensitive to the assumptions regarding the cross-sections of the 
proton capture reactions by the heavy magnesium isotopes and 
by 26 Al: note that the largest enrichment, [Al/Fe] = 1 is con- 
sistent with the aluminum content of the ejecta of the largest 
masses studied here, which are also those mostly depleted in 
oxygen. The production of aluminum is associated to mag- 
nesium burning, which we predict to be sligthly depleted, in 
agreement with the findings of some research groups. 

The sodium yield is the most uncertain, due to the huge 
uncertainties associated to the reactions of the Ne-Na cycle. 
Our models confirm that when the maximum allowed values 
for the reaction rates of the proton capture reaction by 22 Ne nu- 
clei are adopted, the most oxygen poor ejecta are also sodium 
rich, but the exact extent of the sodium enrichment, and the 
confirm that a clear anticorrelation exists, can be hardly fixed 
with the present cross-sections: the theoretical uncertainties re- 
lated to the sodium content amount to 0.6 dex. 

The yield of lithium and fluorine are most sensitive to the 
mass loss rate adopted. With our standard choice, our models 
predict a O-Li trend which is in excellent agreeent with a recent 
investigation based on the lithium content of TO stars in NGC 
6752. The fluorine content is expected to be extremely poor in 
any case, the exact abundance being determined by the details 
of the mass loss description. 
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